Quasiparticle-quasiparticle Scattering in High T c Superconductors 
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The quasiparticle lifetime and the related transport relaxation times are the fundamental quantities 
which must be known in order to obtain a description of the transport properties of the high T c 
superconductors. Studies of these quantities have been undertaken previously for the d-wave, high- 
T c superconductors for the case of temperature-independent elastic impurity scattering. However, 
much less is known about the temperature-dependent inelastic scattering. Here we give a detailed 
description of the characteristics of the temperature-dependent quasiparticle-quasiparticle scattering 
in d-wave superconductors, and find that this process gives a natural explanation of the rapid 
variation with temperature of the electrical transport relaxation rate. 

PACS numbers: 74.25.Fy 



Early measurements of the surface impedance of 
the high temperature superconductor YBa2Cu306+a; 
(YBCO) at GHz frequencies |1[ and at THz frequencies 
M found that the real part of the conductivity, o\ (T) , ex- 
hibited a strong peak as a function of temperature when 
the temperature was lowered below the critical temper- 
ature T c . This effect was interpreted as being due to 
a rapid increase in the transport scattering time of the 
superconducting quasiparticles as the temperature was 
lowered. The rapid increase in the scattering time below 
T c is confirmed by Hall-effect measurements in the flux- 
flow regime 0] and by thermal Hall-effect measurements 
B , and is now well established (further and more recent 
evidence is reviewed in pUfl]). 

Obtaining a quantitative measurement of the tempera- 
ture dependence of this transport relaxation time has not 
been easy, and it is only with the measurements of Hos- 
seini et al. B that information sufficiently precise to test 
current theoretical ideas has become available. These 
recent measurements show that the transport relaxation 
rate is essentially independent of temperature below 20K, 
and increases at least as rapidly as T 4 above this temper- 
ature. In comparing their results with the most relevant 
of the current theories, Hosseini et al. found that their 
T 4 experimental result for the relaxation rate was about 
one power of T faster than the T 3 relaxation rate ob- 
tained in the theory of quasiparticle scattering by spin 
fluctuations in a model for d x 2_ y 2 superconductivity pi. 

The quasiparticle relaxation time is the mean free 
time between collisions of a quasiparticle. The electrical 
(or thermal) transport relaxation time is, roughly, the 
mean free time between those collisions that significantly 
change the electrical (or heat) current. Understanding 
these relaxation times and their differences is central to 
understanding the transport properties of superconduc- 
tors (e.g. see fi). Quasiparticle scattering by impurities 
(relevant at the lowest temperatures) has been studied 
intensively (representative references are [pi and 0) and 




has been found to lead to a number of unusual properties, 
including the phenomenon of "universality" predicted by 
and demonstrated experimentally by Taillefer et 
. Inelastic quasiparticle scattering has been much 
less intensively studied theoretically, and the one relevant 
theoretical study M which does exist does not appear to 
give a sufficiently rapidly varying relaxation rate at low 
temperatures, as noted above. It is therefore of inter- 
est to investigate further different possible mechanisms 
of inelastic quasiparticle scattering. 

In many heavy fermion metals, the electrical resistiv- 
ity p(T) at very low temperatures is found to vary as 
p{T) = p Q + AT 2 . According to [jn) such a T 2 tem- 
perature dependence is usually taken as a criterion for 
the identification of Fermi- liquid behavior |L2] , whereas 
p3| notes that this T 2 dependence could also arise from 
scattering from spin fluctuations. In any case, the fact 
that serious cases have been made that the quasiparti- 
cle lifetime in heavy fermion metals might be limited ei- 
ther by quasiparticle-quasiparticle scattering (the Fermi 
liquid interpretation) or by scattering from spin fluctu- 
ations, suggests that both these mechanisms should be 
investigated in the case of the high T c (d-wave) supercon- 
ductors. Furthermore, ARPES studies of high T c super- 
conductors have been interpreted as giving evidence that 
the ARPES line widths are linked with electron-electron 
interactions |L4| . Because scattering by spin fluctuations 
has already been investigated for d-wave superconductors 
J7J (as well as in s-wave superconductors @]), our article 
is devoted to the study of quasiparticle-quasiparticle scat- 
tering. Interestingly, although both give a T 2 tempera- 
ture dependence in the low temperature limit p^Jl^| for 
a normal metal, we will find that the predicted tempera- 
ture dependences are different for the case of a d-wave su- 
perconductor. These two mechanisms should thus be ex- 
perimentally distinguishable in d-wave superconductors. 

It should be emphasized that the superconducting 
state of the high T c superconductors is by no means 



well understood. In looking at the inelastic scattering 
of quasiparticles in this state (assuming that quasipar- 
ticles exist) it is desirable to get as broad a view as 
possible of a number of different potential mechanisms 
for such scattering (not only the spin fluctuations and 
quasiparticle-quasiparticle scattering just mentioned, but 
also order parameter phase fluctuations Jig , stripe fluc- 
tuations |L7| and phonons) before reaching definitive con- 
clusions. Another point of interest is that even if spin 
fluctuations are an important source of inelastic scat- 
tering at high temperatures, the spin susceptibility is 
expected to decrease in the superconducting state, and 
the quasiparticle-quasiparticle scattering should then be- 
come more important relative to the spin fluctuation 
scattering as the temperature is lowered. The character- 
istics of quasiparticle-quasiparticle scattering elucidated 
in this article, and in particular its agreement with exper- 
iment, suggest that it has considerable promise as an ex- 
planation of the low temperature temperature-dependent 
transport properties. 

According to Hosseini et al. p], the rapid tempera- 
ture dependence of the transport relaxation rate observed 
at low temperature would be expected in any situation 
where the inelastic scattering comes from interactions 
that are gapped below T c . The end result of our low 
temperature calculation (described below) of the trans- 
port relaxation rate for the electrical conductivity result- 
ing from quasiparticle-quasiparticle scattering, namely, 



tIi X = f(T)exp(-Au/k B T) 



(1) 



has just such a gap, making quasiparticle-quasiparticle 
scattering an attractive possible explanation of the low 
temperature inelastic scattering in YBCO. Here, the gap 
Ay is some fraction of the maximum superconducting 
gap, and f(T) is a prefactor which is relatively slowly 
varying, but nevertheless important for the fit of the ex- 
perimental data. 

Quasiparticle-quasiparticle scattering at low temper- 
atures in d-wave superconductors has some interesting 
and novel properties. The scattering of nodal quasipar- 
ticles yields a quasiparticle relaxation rate varying with 
temperature as T 3 . For a typical Fermi surface corre- 
sponding to an optimally doped CuC>2 plane of YBCO 
p8[ , however, such processes are all normal processes 
(conserving the total momentum with no added recip- 
rocal lattice vector) and so do not contribute to the re- 
laxation rate observed in electrical transport (e.g. see 
|fi2|). The processes that determine the electrical trans- 
port relaxation time are the quasiparticle-quasiparticle 
umklapp processes, and these are forbidden unless the 
energy of one of the incoming quasiparticles is greater 
than a threshold energy, Ajj. This is the reason for the 
exponential dependence on Ajj occurring in Eq. nl It will 
be shown below that Eq. FP, which is characterized by the 
exponential factor that varies rapidly with T at low tem- 



peratures, gives good agreement with the experimentally 
measured temperature dependence of rj . 

Calculation of the Quasiparticle Relaxation Rate. The 
quasiparticle lifetime r can be evaluated using the 
"golden rule," and is given by 
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Here M klk2kski is a matrix element that will contain BCS 
coherence factors since we are treating the scattering of 
BCS-like quasiparticles. This matrix element contains a 
5 function conserving the quasimomentum such that 



ki + k 2 = k 3 + k 4 + G 



(3) 



where G is a reciprocal lattice vector. The processes 
for which G = are called normal processes, whereas if 
G^O the processes are called umklapp processes. Also, 
n k = n °(Ek) = {exp(E k /k B T) + 1) _1 is the equilibrium 
value of Fermi Dirac distribution function n k , and E k is 
the quasiparticle energy. 




FIG. 1. The Fermi surface associated with a single C11O2 
plane of YBCO. The superconducting gap varies with mo- 
mentum along the Fermi surface, going to zero at the nodes 
(indicated by solid circles) . The wave vectors Qi on the Fermi 
surface are associated with quasiparticles involved in an umk- 
lapp process (see text) and G is a reciprocal lattice vector. 

At temperatures much less than the maximum gap, the 
thermally excited quasiparticles have momentum vectors 
lying close to the gap nodes on the Fermi surface (see Fig. 
PJ) . In the scattering of a thermally excited quasiparticle 
by another thermally excited quasiparticle the outgoing 
quasiparticles must also have momentum vectors lying 



close to the gap nodes in order to conserve energy. It 
can be seen by studying the Fermi surface geometry of 
Fig. |] that the scattering processes in which only nodal 
quasiparticles are involved must be normal processes. 

Now, the current associated with quasiparticles lying 
close to the nodes in a d-wave superconductor is given by 
the expression 



e rife, 

m* 

k 



J = 2>— » 



(4) 



i.e. the quasiparticle current is proportional to the to- 
tal quasiparticle momentum (e.g. see Q). Because the 
scattering processes just discussed are normal processes, 
they can not change the total quasiparticle momentum, 
and hence can not contribute to the transport relaxation 
rate associated with electrical current. 

It is easily seen that the normal processes discussed 
above cause significant changes to the heat current car- 
ried by the quasiparticles. These normal processes thus 
determine both the quasiparticle relaxation rate and the 
transport relaxation rate appropriate for the quasiparti- 
cle contribution to the thermal conductivity. Therefore 
we comment briefly on their temperature dependence. 
The excitation energies of the nodal quasiparticles can 
be parameterized in the usual way as 



E k = \/(vfPi) 2 + O2P2) 2 



(5) 



where the momentum p is measured from the node and 
has components j>2 along the Fermi surface and j>\ per- 
pendicular to it. The matrix element M in Eq. is taken 
to be independent of momentum (except for the 5 func- 
tion conserving momentum), and Eq. is used. A scal- 
ing argument applied to the momentum integrations then 
yields the result 



DT i 



(6) 



for the temperature dependence of the quasiparticle re- 
laxation rate at temperatures well below the energy gap 
(D is a constant). This same result (except for a change 
of the constant D) is obtained when the appropriate d- 
wave BCS coherence factors are included in the matrix 
element. 

Quasiparticle-quasiparticle umklapp processes do 
change the total quasiparticle momentum and electrical 
current and hence determine the electrical current trans- 
port relaxation time r e i. The quasiparticle momentum 
vectors for one particular umklapp process are shown in 
Fig. fil. Here, the momentum Qi of the low energy quasi- 
particle whose relaxation rate we wish to calculate lies in 
the vicinity of a node. The vector Q 2 is determined by 
the parallelogram construction indicated in Fig. fil and 
by the fact that it lies on the Fermi surface. The four 
quasiparticle momenta satisfy Qi + Q2 = Q3 + Q4 + G 
where G is the non zero reciprocal lattice vector indi- 
cated. A study of Fig. fi] shows that the quasiparticle 



Q2 is the lowest energy quasiparticle that can enter into 
a collision with the nodal quasiparticle Qi in an umk- 
lapp process. The energy of the quasiparticle at Q2 is 
called Ajj (U for umklapp). We expect that the umk- 
lapp process scattering rate will be proportional to the 
mean number of quasiparticles in a state of wave vector 
Q2, which is exp(— Ajj/ksT) for k B T <C Ajj. Umklapp 
processes involving collisions with a quasiparticle with its 
momentum and energy fairly close to those of quasiparti- 
cle Q2 occur for quasiparticles in the neighborhood of Q2 
shown by the shaded region in Fig. fil. The sum over all 
of these umklapp processes gives a relatively slowy vary- 
ing temperature dependent prefactor to the exponential 
temperature dependence just mentioned, as we will now 
indicate. 

For kj in the neighborhood of Qi, let Tiki = hQi + Pi- 
The quasiparticle energy for k2 in the neighborhood of 
Q2 can then be written, in a manner similar to Eq. |g, as 



E k2 = Au + v' 2 p 2 + 



F 2 
Pi 



2A 



(7) 



with a similar equation for E ki . With this parameteri- 
zation the the integrals over the momentum and energy 
conserving 8 functions in Eq. @ can be done analytically, 
giving the result of Eq. Ill Here the prefactor /(T) is (to 
within a constant) 
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dy I d9E(x)F(y)G{x,y,9) (8) 


with the function E(x) = exp(—x 2 j 2 A[j/k B T), F(y) = 
exp[-yAu/(2k B T)}, and G(x,y,0) = (Z - 6)(1 - 
n°{u))/{aZ). Also, a{6) = ( 1 ~ 1 cos9 - sin9) 2 , b(x,ff) = 
a + [(7V7) — x]cos6 + (7' + r yx)sin6, Z — y/b 2 + ay, 
and u = 2Aj/(Z — b)/a. These equations contain four 
undetermined parameters, a — vf/v' f ,j = V2/vf,Y = 
v' 2 /(-/2v'p) and A\j. For an initial investigation of the 
integral for f(T) we chose 7 = 1/14, in agreement with 
experiment fll9[ |, and also made the arbitrary choices 
a = 1,7' = 1/(14V2), and A u /k B = 105K. We find 
that, in the temperature range of interest (20K < T < 
60K), f(T) = CT 2 to an accuracy of about 2%. This 
approximate T 2 temperature dependence is not sensitive 
to reasonable variations of the parameters. Thus, to a 
good approximation 



CT 2 n°(A [/ )[l-n°(A c/ )]. 



(9) 



In this last result, the exponential function has been re- 
placed by the product of Fermi Dirac distribution func- 
tions, which should give a somewhat more accurate result 
as the temperature is raised. With the parameterization 
of Eq. R, however, this result is still correct only in the 
limit k B T < A v . 

The electrical transport relaxation rate r^ 1 has been 
determined experimentally H by fitting the microwave 



conductivity determined at a number of different frequen- 
cies to a Drude line shape (which it fits well). The ex- 
perimentally determined values of t~ { are reproduced in 
Fig. g. For the chosen value 105K of Ajy/Zcs, the the- 
oretical result of Eq. can be seen, in Fig. @, to be in 
agreement with experiment to within the experimental 
error. From the Fermi surface shown in Eq| , from the fact 
that the superconducting gap is found to vary roughly as 
\cosk x — cosk y \ on the Fermi surface J18| , and from the 
parallelogram construction of Fig. m, we find that Ajy is 
about two thirds of the maximum superconducting gap, 
A max . Based on this reasoning our value of Ajj/ks of 
105K yields a A max of 14 meV. The value A max is not 
a very well established value experimentally. For exam- 
ple, for the various samples studied in fll8fl , A max has 
values which lie between approximately 11 meV and 31 
meV (including the uncertainty due to experimental er- 
ror). Given the uncertainties in our method of estimating 
A ma x from our scattering rate formula, and the fact that 
the value that we do obtain is within the bounds estab- 
lished by the ARPES results, the agreement of our theory 
with experiment must be considered satisfactory. 
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2. A comparison of the theoretical curve of 
versus temperature as determined by Eq. PJ with 
0.0139 x 10 11 K~ 2 sec' 1 and Au/k B = 105K, and the 
experimental results (with error bars) of Hosseini et al. |5J . A 
constant has been subtracted from the experimental data so 
that the low temperature limit of tJ 1 is zero. 



d-wave gap, there is an energy threshold for these umk- 
lapp processes such that one of the incoming quasiparti- 
cles must have an energy greater than a threshold energy 
Ay (Ay is some fraction of the maximum superconduct- 
ing gap). This gives the electrical transport relaxation 
rate t~[ a T 2 exp(—A(j/kBT) temperature dependence 
at low temperatures. This theoretical result reproduces 
the rapidly varying temperature dependence of rj 1 ob- 
served H at low temperatures, as can be seen in Fig. ||. 
Quasiparticle-quasiparticle scattering is thus a promis- 
ing mechanism for understanding the low temperature 
inelastic quasiparticle scattering rates in d-wave super- 
conductors. 
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Conclusions This article gives a detailed description of 
the characteristics of quasiparticle-quasiparticle scatter- 
ing in a high T c superconductor such as YBCO, valid at 
temperatures well below the maximum energy gap. The 
quasiparticle relaxation rate and the transport relaxation 
rate appropriate for a description of the microwave elec- 
trical conductivity are found to be controlled by different 
processes. The quasiparticle relaxation rate is due to the 
scattering of nodal quasiparticles off one another and has 
a T 3 temperature dependence. The electrical transport 
relaxation rate on the other hand is due to quasiparticle- 
quasiparticle umklapp processes. In the presence of the 



